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ABSTRACT 


The  present  status  of  the  one-dimensional  computer  routine 
is  discussed  including  recent  modifications  in  the  grain  burning  phase 
of  the  reaction  process.  Calculations  of  minimum  pressure  for  initiation 
of  TNT  are  presented  and  compared  with  experimental  values.  Use  of  such 
comparisons  for  determining  the  rate  parameters  of  the  ignition  process 
are  considered. 

Calculations  are  also  presented  on  the  detonation  behavior  of 
a  typical  composite  double-base  propellant.  Minimum  pressures  for  ini¬ 
tiation,  reaction  and  pressure  profiles  through  the  wave,  and  wave  tra¬ 
jectory,  are  shown.  The  material  is  found  to  be  marginal  for  support 
of  a  steady  state  detonation. 

The  theoretical  backgound  for  the  planned  experiments  on  the 
equation  of  state  of  solids  is  discussed. 
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SECTION  1 

INTRODUCTION 


This  report  presents  the  most  recent  developments  in  the  one¬ 
dimensional  treatment  of  detonation.  In  particular,  the  procedures  for 
handling  the  energy  distribution  between  solid  and  gaseous  phases  now 
appear  to  be  satisfactory.  Provision  has  been  made  for  two  grain  burn¬ 
ing  processes  which  may  conform  either  to  a  surface  ignition  regressive 
grain  burning  model,  or  to  a  point  ignition  progressive  grain  burning 
model.  It  is  next  planned  to  incorporate  these  provisions  into  the 
two-dimensional  computer  routine,  and  when  this  is  accomplished,  it  is 
believed  that  the  procedures  for  pressure-initiated  detonation  will  be 
essentially  complete.  During  the  remainder  of  the  program  it  is  planned 
to  undertake  calculations  designed  to  investigate  the  conformity  between 
computed  and  experimental  behavior.  An  important  part  of  this  latter  step 
is  the  work  on  experimental  study  of  the  equation  of  state  of  propellant 
and  explosive  materials  which  is  presently  in  progress. 

Work  on  thermal  initiation  of  detonation  has  been  held  in 
abeyance  pending  the  development  of  operational  procedures  capable  of 
providing  acceptable  computer  running  times. 
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SECTION  2 

MATHEMATICAL  DEVELOPMENTS 


Two  modifications  were  made  to  the  one-dimensional  computer 
routine  during  this  report  period.  The  first  of  these  is  concerned 
with  the  computation  of  energy  distribution  between  solid  and  gaseous 
phases  in  the  detonation  waves;  the  second  is  concerned  with  adaptation 
to  a  point  ignition  internal  grain  burning  model  which,  as  discussed 
in  the  last  progress  report,  appears  to  be  more  pertinent  to  the  detona¬ 
tion  of  homogeneous  materials  than  the  external  surface  burning  model 
presently  used. 

The  problem  of  energy  distribution  between  phases  was  dis¬ 
cussed  to  some  extent  in  the  last  quarterly  report.  The  specific  methods 
used  to  solve  this  problem  were  found  to  be  unsatisfactory  because  the 
solution  became  unstable  a  short  time  after  ignition  of  the  charge.  This 
instability  was  subjected  to  considerable  study  during  the  reporting 
period,  and  was  found  to  be  caused  by  a  coupling  of  small  oscillations 
in  the  values  of  v  ,  e  and  p.  In  order  to  avoid  this  behavior,  it 
now  appears  necessiry  to  compute  v^,  Vg  by  the  Integration  of  equations 
defining  the  relation  between  dvg ,  dvg  and  dv.  These  are  derived  as 
follows. 


Considering  first  the  dvg,  the  term  is  assumed  to  depend  upon 
the  sum  of  two  effects,  the  first  to  be  derivable  from  the  compression 
or  expansion  of  a  mixture  of  solid  and  gas,  with  simultaneous  transfer 
of  material  from  the  solid  to  the  gas  phase,  the  second  to  be  derivable 
from  the  addition  of  heat  to  the  gas  phase,  at  constant  average  volume. 
(Some  changes  have  been  made  in  the  nomenclature  used  in  this  program. 
New  definitions  are  given  at  the  end  of  the  report.) 
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First,  one  has: 

de  ■  -  pdv  (1) 

s  s 


From  the  solid  equation  of  state: 

(P  +  -  (  7,  -  1)  e,  +  C  (2) 

(P  +  dp  -  (  7^  -  1)  de^  -  -  (  7^  -1)  pdv^ 


Then: 


/3 

(  _—  +  7  )  pdv  -  -  V  dp 

P  8  8  8 


412E - (-^  +  V) 

din  p  s 


And  for  the  gas  O  ■  0) 


iine_  . 

din  v_ 


Also: 


(1  -  f.)  +  f,*.  .  V 


(l-f)dv  -  vdf  +fdv  +vdf  ■dv 

'  s'g  gs  ss  ss 


(3) 


(4) 


(5) 


From  (3)  and  (4) 


s 


(  \  +  ^/P  ) 
7 

g 


Combining  with  (5)  and  solving  for 
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a  m 


dv' 

g 


dv  -  (v  -  V  )  df 

'  s  a'  s 

1  -  f 

s 

1 - . 

V  Y 

1  •  * 

Vg(  7g  +  ^  /P) . 

A  DIVISION  OF 

(6) 


The  effect  of  heet  addition  to  the  gee  phase  at  constant  average  volume 
haa  been  discussed  In  previous  reports.^  Considering  the  following 
diagram  for  one  gram  of  material: 


gas  phase 

- 1 - 

solid  phase 

*vol  ■  V  f 

1 

vol  ■  V  f 

g  g 

1 

s  s 

f  de  ■  dE  -  f  pdv 
g  g  g  g 

f  de  ■  pf  dv 
s  s  s  s 


From  the  equation  of  state: 

^*s  *  7  ~-i  ''s 

de  ■  — ^—r  P<lv  +  V  dpi 

g  7^  -1  [  g  g  J 

Also: 


I  -  f 


dv 


dv 


Elimination  of  dp,  de  ,  de  and  dv  between  equations  7-11  gives 

s  g  s 

(7„  -1)  dE 

dv  "  ■  *  ■  ''  *> 

g  (I  -  f  )  P7»  +  (p7.  +/3)  V, 


(1  -  £,)' 


(7) 

(8) 


(9) 

(10) 


(11) 


(12) 


One  then  has 

dv  ■  dv  '  +  dv 
g  g  g 
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By  similar  methods,  the  expression  for  dv^  Is  found  to  be; 


dv  -  (v  -  V  )  df 

dv  - -S - 1 - S- 

S  V 

f  +  (1  -  f  )  -*- 

8  s'  V 


(  7,  +  i^/P) 


(13) 


and 


de  ■  -pdv 
a  s 


In  order  to  treat  materials  with  a  self -reactive  binder,  an 
additional  grain  burning  equation  has  been  added  to  the  scheme.  Identical 
In  form  with  the  one  presently  used  but  having  different  rate  parameters 
and  associated  with  a  different  heat  of  reaction.  This  requires  an 
additional  term  In  the  calculation  of  the  chemical  heat  release.  The 
previous  differential  equations  are  then  differenced  to  give  the  follow¬ 
ing  computational  procedure: 

IGNITION  PHASE 

Treated  as  before 


GRAIN  BURNING  PHASE 

Set  Initial  conditions : 


eo  •  e° 
s 

f°  -  1  ;  fj  -  1-Fj 


1  ;  fj  -  1-F, 


-  1-F, 


(14) 


(ff  b  1/3 


(fj  ) 


1/3 


B^  At  p 
p  ”  /w  1 


exp 


- - 1 

U:  /  vJ 


(15) 
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1/3 


(£3" 


B3  At  P  f  "  ^3 

f”'*  ♦  1  “"I  +1^'  /-/J 
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(16) 


■  0.8  exp  (-d”^S+0.2  exp  (-Ad""*"^) 


(17) 


_n+l_._  Do  V  _2  /..n+l  . 

D  "(D  +  “  )  TT  (t  -  t) 

“  P 


(18) 


-n+1  n+l  n+1 

■  '2  '2  +  h  ‘3  *  ^ 


(19) 


n+1  -n 


(V  V  ■  )  -  (v"  -  v"  X^*  - 

g  »  B  »  ' 


n+1  n 

V  ■  V  +  ^ , 

8  8  ^n+1 


n+1  n 


V  + 

8 


*  <1  - 

n+1  -  n.  .  n  n  .  .,n  ,n+l  . 

(v  -v)-(v  -v)(f  -f  ) 

J _ '  g  _*  *  -  * 

1-f.fi-  :i-  Jm _ ^1 

I  7"  +^/p"J 


(20) 


(21) 


(y;  -1) 

(e;  -.g")(f;  +  F2Q2(f2" 

■ 

)+  F3Q3 

,(£3"  -f3"^S40^(f^" 

(1- 

£”♦*) 

8 

P"7j+  (p"7"  +^) 

V 

V 

8 

n  /I  -  f^"-^l  1 

fn+l 

n+1 


n+1 


b 

n+1 


7  +  — 

.0 

V*  j^,n+l 


» 

b  \2 


n+1 


(lo-  1) 


(22) 


(23) 
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.1) 


[fO+l 

■<=>  +  ■ 

A],t- 


«2  )+  ‘^3  > 


.  _  ,n+l  ,  ^  n  .  «  n+l/2v  »  v  v 

+  %  ’h  )  ■  (P  +  2q  '  )  <  - 2 - 


1  g  ‘  s 

n+1  n+1  .  n+1  .  v  -  v  ''.I 

f  V  +  (  7„  -  1)  (  ) 

- 1  1 

1  s  s  g  2  I 

(24) 


-  n+1  -  n  .  n  .  n+1  .  «  n+l/2v  ,  v  v  ^  ,,n  ,n+l  v 

e  •  e  -  (p  +  p  +  2q  )  (  - 5 -  +  F2Q2  (fj  "  ^2  ^ 


—  ^  .n+l  .  _  . ,n  _n+l. 

+  F3Q3  (£3  -  £3  )  +  (£^  -  £^  ) 


(25) 


n+1  n  .  ,  n  .  n+1  .  .  n+1/2.  .  *  v”  ^ 

**  •«,+(?+?  +  2q  '  )  (  _1 _ S_) 


(26) 


n+1  n+1 

P 


(27) 


Other  details  o£  the  program  remain  the  same. 

The  second  modl£lcatlon  concerning  adaptation  o£  the  routine 
to  an  Internal  grain  burning  model  has  turned  out  to  be  relatively 
stralght£orward.  It  was  shown  In  the  previous  progress  report^  that  the 
grain  burning  rate  laws  £or  the  two  models  are  Identical  If  £  Is  defined 
in  the  external  burning  case  as  the  fraction  of  material  unburned,  and 
In  the  Internal  burning  case  as  the  fraction  of  material  burned.  It 
follows  very  readily  that  conversion  from  the  former  to  the  latter  problem 
can  be  accomplished  by  changing  the  sign  of  Q2>  Q3(  82  and  83  In  the  Input, 
and  redefining  the  Initial  values  of  the  f's  as  follows; 


f2°-0 


f3°.0 


'I-*'! 


(28) 
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Integration  la  terminated  when  £  >  1. 

The  proper  procedure  for  la  aomewhat  In  doubt  with  thla  model. 
However,  the  dlffualon  rate  law  la  algnlf leant  here  only  In  the  caae 
of  cooq^oalte  double-baae  propellanta  with  grain  alae  of  the  order  of 
5  mlcrona  or  leaa.  If  It  can  be  aaaumed  that  the  centera  of  Ignition 
are  located  In  the  Interface  between  oxldlaer  grain  and  binder,  then 
the  moat  reaaonable  procedure  la  to  uae  the  aame  geometry  for  the  dlf¬ 
fualon  proceaa  In  both  modela.  Although  thla  may  falalfy  to  aome  extent 
the  form  of  the  dlffualon  reaction  profile  through  a  detonation  wave. 

It  would  not  algnlflcantly  affect  the  reaction  zone  length  and  ao  should 
not  appreciably  alter  detonation  behavior. 
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The  one -dlnene tonal  program  In  the  form  deecrlbed  in  Section  2 
hae  been  used  to  examine  the  detonation  behavior  of  TNT  and  of  a  repre¬ 
sentative  composite  double-base  propellant.  Input  data  for  the  two 
materials  Is  as  shown  In  Table  I  except  as  noted  In  the  Illustrations. 

The  external -regressive  grain  burning  formulation  of  the  grain  burning 
reaction  was  used  In  these  cases. 

Of  particular  Importance  In  calculating  the  detonation  behavior 
of  real  materials  Is  the  selection  of  proper  values  for  the  rate  parameters 
A]^,  A2i  and  82-  Data  for  TNT  are  available  from  two  sources.  The 
first  of  these  consists  of  thensal  decomposition  studies;  this  source  has 
been  reviewed  and  used  by  Cook.^*^  The  second  consists  of  experimental 
work  carried  out  at  the  Naval  Ordnance  Laboratory  on  the  shock  pressures 
necessary  for  Initiation  of  detonation  In  TNT.^  Such  pressures  are 
strongly  dependent  upon  the  values  used  for  Bi  and  Ai  and  the  data  can 
therefore  be  used  to  evaluate  these  parameters  If  the  other  needed  Infor¬ 
mation  about  the  material  Is  available.  As  has  been  discussed  In  previous 
reports,  the  assumption  that  the  Ignition  zone  Is  a  high- temperature  region 
In  the  vicinity  of  a  bubble  or  void  which  has  been  passed  over  by  a  shock, 
leads  to  the  postulate 

A^  -  k^n 
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TABLE  I 

INPUT  DATA  FOR  CALCULATION 


TNT 

Common 

Composite  Double  Base 

C  -  795  cal/ga 

a  -  4.0  X  10’® 

C  ■■  664  cal/gm 

/3  -  SO  kbar 

7  -  1.3 
'0 

/3  ■  SO  kbar 

-  1130  cal/gm 

b  -  1.03 

■■  622  cal/gm 

Qj  ■  Q2  ■  1130  cal/gm 

-4  2 

■  10  cm  /sac 

Qj  ■  Q2  ■  380  cal/gm 

Q3  -  Q3  -  0 

D  ■  10  ^  cm^/sec 

0 

Qj  -  Qj  ■  864  cal/gm 

Q4  -  0 

w  >  20  kbar 

Q4  -  0 

Aj/R  -  4715®K 

m  •  1 

Aj/R  -  5500®K 

Aj/R  -  9430°K 

Aj/R  -  11,000°K 

-  1.08  X  10^^ 

Aj/R  -  19,400®K 

Bj  -  1.08  X  10^^ 

Bj  -  1.6  X  10^^ 

Bj  -  0 

12 

Bj  -  1.6  X  10 

•  0.01 

Bj  -  4.6  X  10^^ 

Fj  -  0.99 

Fj  -  0.01 

0 

i 

Fj  -  0.205 

F^  -  0 

F^  -  0.470 

Z  ■  0.14  cm 

Fj  -  0.315 

■  1.5  gm/cm^ 

Z^  ■  2.5  X  10  ^  cm 

...-3 

P  -  93.8  kbar 

Zj  -  12  X  10  cm 

tp  ••  1.5  ^sec 

2 

-  1.80  gm/cm 

P  -  100  kbar 

t  p  ■  4  ^  sec 
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The  procedure  for  obteining  velues  of  the  reCe  peremeters  le 
to  compere  computed  velues  of  the  minimum  pressure  for  Initietion  of 
detonetlon  obteined  with  different  velues  of  end  with  the  experi¬ 
ments!  velue.  The  computed  velues  ere  obteined  by  setting  B2  end  B3 
equel  to  zero  so  thet  the  greln  burning  reectlons  ere  suppressed,  end 
by  using  e  low  velue  of  Fx  so  thet  the  Ignition  reection  elone  Is  not 
cepeble  of  supporting  e  detonetlon  weve.  The  celculetlon  Is  Inltleted 
by  en  Imposed  shock  of 'v  100  kbers  heving  e  duretion  of  1.5  fxsec,  which 
is  heevy  enough  to  ceuse  Ignition.  This  dlsturbence  then  propegetes 
Into  the  meter lei  es  e  fedlng  weve,  end  the  meximum  weve  pressure  et  the 
lest  zone  which  Ignites  Is  the  minimum  pressure  for  Inltletlon. 

Figure  1  shows  the  trejectory  of  such  e  fedlng  weve  obteined 
trlth  velues  of  end  B^  derived  from  the  dete  of  Cook.  This  gives 
directly 

A2  ~  43,400  cel/mol 
B2  •  1.08  X  10^^ 

whereupon 

Aj  -  21,700  cel/mol 
B^  -  1.08  X  10^^ 


For  purposes  of  the  celculetlon  B2  end  B3  wre  set  equel  to  zero.  Dete 
on  the  heet  of  explosion,  density,  end  specific  best  of  TNT  were  obteined 
from  Cook.^  Velues  of  the  other  input  peremeters  were  epproxlmetlons 
obteined  es  discussed  In  previous  reports. 

Comperlson  of  the  fede  point  with  the  curve  of  weve  pressure 
shows  thet  fedlng  occurs  et  e  pressure  of  <v90  kbers.  This  velue  Is 
somewhet  uncerteln  beceuse  of  the  rether  wide  fluctuetlons  which  occur 
in  the  computed  mexlisum  pressure.  Accurecy  could  be  Increesed  to  eny 
desired  velue  by  Increeslng  the  steblllty  conditions  imposed  on  the 
integretlop  end  decreesing  the  mesh  size.  This  would,  however,  incresse 
the  running  time,  end  wes  not  deemed  to  be  of  importence  to  the  present 
celculetions .  The  experlmentel  velue  reported  by  NOL  Is  ebout  30  kber. 
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FIGURE  1.  POSITION  AND  MAGNITUDE  OF  PEAK  PRESSURE  IN  FADING 
WAVE  IN  TNT  (A,/R  -  10,700®K  -  A-ZR  -  21,400®K) 
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It  is  deduced  that  either  the  rate  parameters  obtained  from  thermal  decompo¬ 
sition  data  predict  reaction  rates  much  lower  than  are  actually  obtained 
in  a  detonation  process,  or  that  the  assumed  equation  of  state  parameters 
lead  to  a  value  of  energy  density  behind  the  lead  shock  much  lower  than 
that  actually  obtained.  At  the  present  time  it  is  not  possible  to  choose 
between  these  two  explanation. 

If  the  value  taken  for  is  9500,  then  the  results  are  as  shown 
in  Figure  2.  The  minimum  shock  pressure  for  initiation  here  agrees  very 
well  with  the  experimental  value.  The  A^  corresponds  to  an  A2  for  TNT  of 
'V  20  kcal  which  is  not  regarded  as  unreasonable. 

Figures  3,  4  and  5  show  additional  results,  using  these  values 
of  A.  and  A..  Figure  3  shows  the  reaction  profile  through  the  wave  at  a 
time^of  12  /isec  after  initiation  when  the  wave  has  propagated  about  5  cm 
into  the  charge.  Complete  reaction  behind  the  wave  is  not  being  achieved 
at  this  point,  and  therefore  it  would  not  be  expected  that  the  steady 
state  velocity  would  have  been  reached.  This  is  shown  by  a  plot  of  wave 
trajectory  in  Figure  4,  which  indicates  a  velocity  of  about  4  tm/fisec. 

It  is  expected  that  upon  further  propagation,  the  wave  velocity  would  in¬ 
crease  and  the  reaction  zone  length  would  decrease.  However,  since  the 
overall  length  was  set  at  10  cm,  it  was  regarded  as  unlikely  that  growth 
to  steady  state  would  occur  before  reaching  the  end  of  the  charge.  This 
wave  will  be  examined  further  by  rerunning  the  problem  with  a  longer 
charge  length. 

Figure  5  shows  plots  of  v  ,  v^  and  e  through  the  wave.  The 
reduced  amount  of  scatter  shown  in  the  points  is  a  consequence  of  the 
imposition  of  increased  stability  criteria  on  the  integration.  These 
results  together  with  those  discussed  above  indicate  that  the  mathematical 
procedures  now  being  used  are  quite  well  behaved. 

Calculation  has  also  been  carried  out  on  the  propagation  of  a 
detonation  wave  through  a  representative  composite  double-base  propellant 
using  the  input  data  indicated  in  Table  !•  obtained  from  information 
provided  by  the  Hercules  Powder  Company.  In  this  problem,  two  grain 
burning  reactions  were  used,  one  for  the  ammonium  perchlorate  oxidizer, 
end  another  for  the  double-base  binder.  As  in  the  case  of  TNT,  the  rate 
parameters  for  the  binder  were  obtained  from  thermal  decomposition  and 
adiabatic  self-heating  data.  Those  for  ammonium  perchlorate  were  obtained 
from  data  published  by  the  Aerojet  General  Corporation  on  the  regression 
rates  of  APC  surfaces  in  their  hot  plate  experiments,  which  have  been  dis¬ 
cussed  in  a  previous  progress  report.^  The  aluminum  was  treated  as  inert. 
Grain  sire  taken  for  the  ammonium  perchlorate  was  the  average  size  of  the 
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FIGURE  2 ,  POSmCMI  AND  MAOIITUDE  OF  PEAK  HIESSURE  IN  FADING  WAVE  IN 
TNT  (Aj/R  -  4715°K  -  Aj/R  -  9430®  K) 
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FIGURE  A.  INITIAL  TRAJECTORY  OF  ISTONATKNI  WAVE  IN  TNT 
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FIGURE  5 .  SOLID-PHASE  ENERGY  DENSITY  AND  SOLID  AND  GAS  SPECIFIC 
VOLIMES  THROUGH  A  DETONATKm  HAVE  IN  TNT 
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material  used  In  the  original  mix;  that  for  the  binder  was  1/2  the  average 
distance  between  APC  grains.  The  rate  parameters  for  the  ignition  reaction 
were  taken  as  those  of  the  ammonium  perchlorate  for  two  reasons.  First, 
under  conditions  of  low  shock  intensity  the  APC  rate  law  gives  a  faster 
reaction  than  that  for  the  binder.  This  rate  would  therefore  be  most  im¬ 
portant  in  the  ignition  process.  Second,  NOL  data  from  card  gap  experi¬ 
ments^  support  this  conclusion.  The  heat  of  the  ignition  reaction 
was  taken  as  the  average  of  the  heats  of  reaction  of  the  APC  and  binder 
since  they  are  assumed  to  contribute  equally  to  the  ignition  reaction. 

Figure  6  shows  the  trajectory  and  peak  pressure  in  a  computed 
wave  initiated  with  a  100  kbar  shock.  Initial  wave  velocity  is  found  to 
be  about  5  mm/pisec  slowing  to  about  4  nm/jLisec  after  10 /x sec.  Peak  pres¬ 
sures  decay  rapidly  from  the  initial  Imposed  value  of  100  kbar  to  about 
50  kbar.  It  is  again  pointed  out  that  the  scatter  shown  in  the  computed 
points  is  an  inherent  consequence  of  numerical  integration  methods, 
which,  for  these  problems,  is  particularly  severe  in  the  region  of  peak 
pressure.  It  could  be  reduced,  as  desired,  by  tightening  the  stability 
criteria  for  the  integration;  however,  at  an  Increase  cost  in  computing 
time. 


Figure  6  indicates  that  this  wave  approached  a  steady-state 
condition.  However,  at  the  point  at  which  the  calculations  were  terminated, 
the  pressure  and  velocity  still  seemed  to  be  decreasing  slightly,  and  it 
is  possible  that  fading  would  ultimately  have  occurred.  The  conclusion 
is  that  the  material  is  marginal  for  detonability.  It  is  of  considerable 
interest  here  that  experimental  data  on  similar  propellant  materials  ob¬ 
tained  at  the  Naval  Ordnance  Laboratory  are  consistent  with  the  conclusions 
to  be  drawn  from  these  calculations  that  a  pressure  of  50  kbar  is  somewhat 
above  the  maximum  pressure  required  for  initiation.^ It  would  be  pre¬ 
mature  to  speculate  as  to  the  significance  of  this  agreement  at  the  present 
time.  However,  an  important  difference  between  these  results  and  those 
discussed  previously  for  TNT  lies  in  the  fact  that  the  rate  data  used  for 
the  ignition  reaction  in  the  present  case  were  obtained  by  the  hot  plate 
technique,  while  those  for  TNT  were  derived  from  thermal  decomposition 
data.  Further  study  of  this  point  appears  desirable  because  it  may  shed 
some  light  on  the  question  of  which  experimental  methods  are  most  suited 
for  obtaining  kinetic  information  appropriate  for  predicting  detonation 
behavior . 


Figure  7  shows  some  data  on  reaction  and  pressure  profiles 
through  the  detonation  wave  in  the  above  composite  double-base  propellant. 
It  is  of  interest  here  that  in  the  surface  layers,  the  binder  has  burned 
completely  and  the  APC  only  to  a  slight  extent.  In  the  interior  of  the 
charge,  neither  binder  of  APC  has  reacted  significantly. 
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FIGURE  6.  TRAJECTORY  AND  PEAK  PRESSURES  OF  DETOHATION  WAVE 
IN  CGHPOSITE  DOUBLE -BASE  PROPELLANT 
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The  explanation  for  this  behavior  appears  to  reside  in  the  high 
activation  energy,  (A2) ,  and  thus  the  high  temperature  coefficient  for  the 
binder  reaction.  Reaction  in  the  surface  is  rapid  because  of  the  high 
energy  deposited  in  the  material  by  the  initial  shock.  However,  the  energy 
released  by  binder  reaction  is  not  sufficient  to  drive  a  shock  of  this  in¬ 
tensity  so  that  following  release  of  the  surface  pressure  the  wave  decays. 

As  a  consequence  of  the  decreasing  pressure  in  the  wave,  the  energy  in  the 
solid  drops  to  a  value  at  which  the  binder  no  longer  reacts  at  a  significant 
rate.  Similar  behavior  is  not  observed  on  the  part  of  the  ammonium  per¬ 
chlorate  reaction  because  its  lotrer  activation  energy  makes  it  less  sensitive 
to  changes  in  solid  energy  density.  Following  the  point  at  which  the  grain 
burning  reactions  subside,  the  wave  is  supported  only  by  the  Ignition 
reaction,  and  this  constitutes  the  observed  steady  state  p. ocess. 
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SBCnON  4 

EQUAXION  OF  STATE  STUDIES 


The  equation  of  state  presently  used  for  the  solid  phase  Is 
similar  In  form  to  that  derivable  by  representation  of  the  solid  as  a 
collection  of  harmonic  oscillators.  If  the  possible  frequencies  are  Vg, 
the  equation  obtained  from  such  a  model  Is  the  following:^ 


(29) 


where 


Sin  V 


:k 


I  din  V 


(30) 


Here  the  summed  terms  represent  the  Internal  thermal  energy  of  the  material 
and  0  Is  a  potential  function  such  that  the  Internal  potential  energy 
Is  given  by 


V 

%  ‘  f  ^ 

CO 

For  present  purposes,  Internal  thermal  energy  Is  meant  to  Include  the 
kinetic  and  potential  energy  associated  with  atomic  end  molecular  vibra¬ 
tion;  and  Internal  potential  energy  represents  that  energy  due  to  atomic 
and  molecular  distortion. 
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At  high  temperature  the  summation  in  equation  (29)  reduces  to 


RT 


z 


e 


(32) 


iihere  |  Is  the  aero  point  energy.  If  all  of  the  7g  are  equal,  substitu¬ 
tion  of  (32)  In  (29)  gives: 


(33) 


Ihe  oscillator  model  upon  which  equation  (33)  Is  based  Is 
closely  approximated  only  by  a  perfect  single  crystal.  However,  the  form 
of  the  equation  la  a  consequence  of  the  condition  that  the  potential  and 
thermal  forms  of  Internal  energy  are  respectively  dependent  upon  the 
volume  only  and  upon  the  tea|>erature  only,  and  are  therefore  separable 
Into  different  terms.  If  It  Is  assumed  that  however  poorly  the  simple 
oscillator  model  may  represent  reality  In  other  respecta,  this  separation 
of  thermal  and  potential  energy  is  valid  for  real  Mterials,  then  the  ^ 
equation  Is  of  general  application,  provided  only  that  the  functions  r 
and  which  now  must  be  regarded  as  somewhat  empirical  In  nature, 
known. 


In  treating  detonation  behavior.  It  Is  more  convenient  to 
express  the  equation  of  state  In  terms  of  total  energy,  since  the  conserva¬ 
tion  equation  refers  to  total  energy.  Adding  equation  (31)  to  both  sides 
of  (33),  one  has  ultimately 


•  V  ^ 

fa*. 

J  37'*’ 


dv  I  V  -  e  + 


m  « 

3  N  k  T  +  ^ 
o 

m  R 


(34) 


Thlsla  nowtheform  of  the  solid  equation  of  state  presently  In  use  where 

V 


A 


y  j  ^ 


(35) 


7g  -  7-1 
C  -7q  [3NkT^  +  e] 
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3MkT^  is  the  thermal  energy  content  at  ambient  tei^perature . 

e  >  total  Internal  energy  referred  to  zero  at  ambient  conditions. 

An  objective  of  the  present  program  is  to  determine  the  functions  and 
7g  the  materials  and  the  condition  of  significance  to  detonation 
phenomena.  These  can  be  obtained  from  p-v-T  data.  Partial  differentia¬ 
tion  of  equation  (33)  gives: 


3Nk  (It  )v 

(36) 

or 

^  |lf|v 

(37) 

y\  1 

Thus 

at  any  v  is  given  by  the  slope  of  the  p-T  curve. 

Substitution 

in  (33) 

then  provides 

The  necessary  p-v-T  data  will  be  obtained  by  static  pressure 
methods,  making  use  of  a  piston-cylinder  device  designed  by  Kennedy^^ 
capable  of  an  ultimate  pressure  in  the  order  of  80  kllobars.  This 
apparatus  has  been  delivered  and  is  in  the  process  of  being  set  up  and 
checked  out.  Some  preliminary  development  of  methods  end  techniques  will 
be  necessary.  It  is  estimated  that  initial  equation  of  state  data  will 
be  available  near  the  end  of  the  next  report  period. 
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NOMENCLATURE 


t  - 

p  ■  pressure 
q  ■  von  Neumenn  "q" 

V  ■>  volume 

u  ■  meterial  velocity 
e  ■  specific  energy 
h  ■  Planck's  constant 
k  ■  Boltsman  constant 
r  >  Ignition  delay 

V  ■  empirical  pressure  constant 
^  ■  equation  of  state  constant 
C  ■  equation  of  state  constant 
7  ■  equation  of  state  constant 
T  ■  temperature 

R  ■  gas  constant 

V  m  vibration  frequency 


0  ■  potential  function 
f  ■  fraction  of  material  unreacted 
F  ■  fraction  of  material  In  tone 
A  >  activation  energy 
B  ■  pre-exponential  constant 
m  >  reaction  order 
Q  ■  heat  of  reaction 
B  ■  chemical  energy  released 
(  >  aero  point  energy 
Cy  ■  specific  heat 
L  ■  charge  length 
Z  •  grain  radius 
p  •  density 
a  ■  molecular  diameter 
P  ■  applied  pressure 
U  •  wave  velocity 


Subscripts 

o  ambient  condition 

1  Ignition  process 

2  oxldlser  process 

3  binder  process 

4  diffusion  process 
g  gas  phase 

s  solid  reaction  phase 
1  enert  phase 
G  Grunelsen  constant 
a  Initial  conditions 
z  final  conditions 

Superscripts 


n  time  mesh  number 
J  space  mesh  number 
o  Ignition  point 
bar  average  value 
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